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I. INTRODUCTION 

The idea of a spherical detector of gravitational waves has been considered almost since the early seventies [Q. 
Very accurate calculations have been performed j|] m order to use a self-gravitating object as a possible detector. 
Basing upon these results other authors || have solved the inverse problem for a free sphere in interaction with a 
gravitational wave. In doing so they have shown an omnidirectional sensitivity and the possibility of testing the 
presence of a longitudinal term in the gravitational signal. However a free spherical resonator is not a good device. 
In fact a set of secondary mechanical oscillators (transducers) is needed in order to convert the mechanical signal into 
an electromagnetic one (see ||, Thorne, p. 330]). In the spherical geometry almost 5 transducers are needed to obtain 
all the informations contained in the motion of its 5 quadrupolar modes. The problem of the interaction of a sphere 
with these secondary resonators has been solved in some particular configuration only in the nineties In the 

configurations envisaged by these authors the omnidirectional sensitivity and the possibility of total reconstruction of 
the signal has been shown, in this way realizing a complete observatory with only one device (cfr |8|-[l0|). However, in 
the calculations performed until now, not enough attention has been paid to the possibility of vetoes against spurious 
events (see for instance JTT|] ) and to the influence of thermal noise in the reconstruction of the signal. In this paper we 
have focus our attention upon these two problems, solving them explicitly for the dodecahedron arrangement of the 
transducers first described by [||. Moreover a Lagrangian which describes the interaction of a sphere with N radial 
transducers on the surface has been written. This is important in order to write the most general equations involving 
the motion of the modes of the sphere and of the secondary oscillators. In this way it will be possible to undertake 
a systematic study of other transducers's configurations, enabling spherical detectors to improve their capabilities in 
the framework of gravitational astronomy Jll],[l2| . 

II. LAGRANGIAN OF A FREE ELASTIC BODY INTERACTING WITH A GRAVITATIONAL WAVE 



An elastic body in interaction with a gravitational wave could be described by ordinary elastic theory if its dimen- 
sions are very small compared to the gravitational wavelength. Gravitation acts as an external field of density force 
(e.g. |) 

P = pR l = pE l3 {t)x 3 (2.1) 

where Eij are the components of the "electric" part of the Riemann tensor Riojo- If u(x, t) represents the displacement 
from the equilibrium position in the point x then (e.g. p3| ) 

p-Qp=vV 2 u+(\ + f i)V(V-u)+ P R (2.2) 
Using the standard approach the solution of ( |2.2| ) could be decomposed as 

u(x,t)=J2Mt)ip n (x) (2.3) 

n 

in which ip n (x) are the eigenfunctions with eigenvalue Lo n of 



f 



H V 2 ^ n (x) + V(V • ip n {x)) = -uj 2 Mx) (2.4) 
P P 

The eigenfunctions ip n (x) are a complete set of orthonormal functions with the scalar product defined as (* means 
complex conjugation) 

< g(r, 6, 0)|/(r, 0, 0) >= JL f g * . f pd \ (2.5) 
where M is the mass of the elastic body. Therefore 

77 / V£ • tPn'Pd 3 * = $nn> (2.6) 



Time dependent coefficients A n (i) satisfies the following forced harmonic oscillator equations 

A n +ij 2 n A n =R n (2.7) 

where 

Rn(t) = ^J C(^) ' R(x,t)pd 3 x (2.8) 

In this way a field equation has been transformed in a set of equations for the coefficients of a function series. Now 
we can write the Lagrangian from which these equations are derived. This Lagrangian should be a real expression of 
the generalized coordinates A n (t) and A*(t). As far as kinetic energy of the elastic body (EB) is concerned, one has 



T EB = \I d 3 X P U 2 = UdY^AnK 



(2.9) 



where last equality follows from ( |2.6| ). Potential energy of both the elastic body and the external forces are immediately 
written as 



Vbb = \MY,<fiA«K (2-10) 

n 

Vext = ~M Y, \AnK + A* n Rn] (2.11) 



III. LAGRANGIAN OF THE FREE SPHERE 

If the elastic body is a sphere with radius R, then (see 0) 

u{r,9,cp,t)=Y J A nm im nml (r,e,4>) (3.1) 

nml 

where 

^ nlm (r,d,cj>) = [a nl (r)n + b nl (r)RV]Y lm (6,<t>) (3.2) 

p dji (q nt r) 1(1 + 1) . 

a n im (r) = Cnim R j h a„ irn R ji (k n ir) (3.3) 

ar r 

bnlm (r) = Cnlmjl (q n ir) + d„i m 4; ( r ji (Kir)) (3.4) 

and A nlm (t) = (-l) m A* nl _ m (t). 

Now we briefly sketch the procedure used by || to take into account direction of the gravitational signal. In the 
laboratory frame O, with origin in the centre of the sphere, the gravitational wave is coming by a given direction. Let 
us consider a reference frame O' rotated with respect to O in such a way that z' axis coincides with the direction of 



2 



propagation. All the primed quantities refer to this last system. The radial component of the surface displacement 
in the primed frame is then given by (see eq. (7)]) 

n' ■ = a nl {R)A' nlm ,{t)Y lm ,{e',4>') (3.5) 

nlm' 

We are interested in the radial displacement in the laboratory frame. Let us suppose that O is rotated with respect to 
O' by the hour angle H (the angle between the projection of z' axis on the plane xy and x axis) and the declination 
5 (the angle between z' axis and the plane xy, diminished by 2 it if obtuse) B. In O we have 



nlm \ m' / 



(3.6) 



where ^m m '(H,S) are the rotation matrix for spherical harmonics (see e.g. @). Let us define 



F nlm =Y,A' nlm ,V^ im ,{H,5) (3.7) 



which are the coefficients of decomposition (3.1) in O. These coefficients generalizes in a natural way coefficients F m 

as defined in || eq. (9)], for every order in indices n and I. More precisely one has F m = F 12m . V^ m ,(H,S) are 
unitary matrices, therefore it is possible to invert eq. (|3.7|) and find A' nlm , as linear combinations of F n i m 



Kim' =J2 V ^m*(H,S)F nlm (3.8) 

III 

The kinetic ene rgy of the sphere (S) as a function of the generalized coordinates in the laboratory frame will be 
given by (see ( |2.9|) ) 

Ts = \m ]T Kim'Klm' = hi F nlm F* nlm (3.9) 

nlm 1 nlm 

As far as the elastic potential energy of the sphere and the potential energy of the external gravitational force are 
concerned one has 

V S = \M "llKlm'Klm' = \M "ll F nl m F* lm (3.10) 
nlm 1 nlm 

Vext = ~2^1 [A'nlm'Rnlm' + A'*i m ,R n lm'} = — 2 [G*nlmFnlm + GnlmF* lm ] (3-11) 
nlm' nlm 

where we have set 

Gnlm = y^^mm'Rnlm' (3-12) 

m' 

which are the coefficients of the decomposition of the gravitational force in O. 



IV. COUPLING WITH TRANSDUCERS 



A. General Case 



Let us consider a set of TV harmonic mechanical oscillators with mass and resonant frequency Ui, placed on 
the surface of the sphere in (R,9i,4>i). Let qi(t) be the generalized coordinates which describe radial motion of each 
oscillator with respect to the surface of the sphere. We make the assumption that mj are small with respect to 
M, the mass of the sphere. In this way we suppose that the presence of secondary oscillators does not modify the 
eigenfunctions of the sphere (in the same way one suppose that the transducer of a resonant bar does not modify the 
eigenfunction of the cylinder). Purpose of this section is to write that part of Lagrangian relative to the motion of 
the resonators and to their interaction with the surface of the sphere. 
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Concerning the potential energy of the transducers (t), the choice of generalized coordinates brings to the expression 



Vi 



It, 



m,iW?qf 



(4.1) 



Moreover, if we call yi the coordinates of the i— th oscillator in the laboratory frame, then kinetic energy of transducers 
is written as 



T t = 



(4.2) 



The relation between inertial and generalized coordinates is given by 



qi =y i -n- u{R,0i,(j>i,t) 



(3.6M0) 



It — 2^ a nl{R)FnlmY 1 



nlra 



For convenience purposes, we define, for every I, the matrix whose components are 



P 



(') 



Y, 



(4.3) 



(4.4) 



The connection between these matrices and the matrix B of ref. S is given by B = P^* With these positions the 
kinetic energy of the transducers is 



. TP * 

mi nlm 



film 



E a niWa nn ,(jV)pWv^>r nh! ,r,: rill . 

nlm n'l'm' 



(4.5) 



In conclusion the Lagrangian of an elastic sphere with TV mechanical oscillators on its surface, undergoing an 
external force is 



L = T S +T t -Vs-V t -V ec 



(4.6) 



where terms involved have been defined in (3.9), ( |4.5| ), (3.10), ( f4.1| ) and ( 3.1 1| ) respectively. Generalized coordinates 
are F nhn , FJ m and q if where I = {1, 2, ...}, m = {-I, -I + 1, I}, and i = {1, 2, .., N}. Since F ni _ m = (-l) m FJ m 
it is sufficient to find Lagrangian equations from qi and, say, F^. Therefore, to simplify the problem, we can take 
suitable linear combinations of the coefficients F n i m in such a way to obtain real quantities. In particular, for every I it 
exists a unitary operator which transforms irreducible tensors of rank I in real vectors having 21 + 1 components. 
We will use lower case letters to indicate those real vectors corresponding to irreducible tensors. One can therefore 
write 



(4.7) 



where /„ ; is a vector whose components are f n i a (a G {1, 2, ... ,21 + 1}), while F n i m (m G {—I, are the 

components of F n i. From now on indices a, b, c run from 1 to 21 + 1. In this way, for instance, components of U™' 
will be given by llim- A possible choice for such unitary operators is the following 



in 



(4.8) 



U, 



1 



2m'+lm — V*rn' m + ( 1) $m' —m) 



where m! G {1, 2, . . . , I}. Performing these unitary transformation Lagrangian is obtained as a function of the real and 
independent generalized coordinates f n i a and qi. With this new set of variables, kinetic and elastic potential energy 
of the sphere, kinetic and potential energy of the transducers and potential energy of the external gravitational force 
acting on the sphere become respectively 
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(4.11) 



Ts = \MY J flia (4-9) 

nla 

Vs = lMj2"llfnla (4-10) 

nla 

i L nio 

+ anl{R)a n n'{R)p { aJp^Unlafn'l'a 
nla n'l'a' 

V^l^m^q* (4-12) 
T4 xt = -M^g„ io /„ ;a . (4.13) 

where the matrix = (p^) is given by 

Obviously the potential energy of the transducers is left unchanged. Lagrangian equations 

d ( dL \ dL 

dt KdfrOa) dfnla ~ 

d_ fdL\ _d£ 
dt \dqi J % 

give therefore the following equations of motion for the modes of a forced elastic sphere in interaction with N harmonic 
mechanical oscillators 

Mf n i a + Mu) 2 nl f rda - ^2 a„i(R)m I ujfp^ ) i q I = Mg nla 

i 

(4.15) 

m l q l +171,^2 a nl(r)Pai fnla + m,ujfq l = 
nla 

B. Quadrupolar Interaction 

In which follows we make the simplifying assumptions that secondary oscillators are identical, their common fre- 
quency is equal to u>\2 (the frequency of the mode n — 1, I = 2 of the sphere), their mass is m and the external force 
acting on the sphere has only I = 2 components. Moreover, for high quality factors (t„; to is the relaxation time of the 
nlm mode) 

I I 1 r , ii^l 

\u>12 ~ u>nl\ » ~ t = min{Ti 2m ,T n i m >\; ^ ^ 2 

therefore if the motion of the transducers is studied in the frequency domain near the frequency W12, only modes with 
n = 1 and I = 2 must be taken into account. To simplify notations we set 

P ai =Pai fl2a = fa 9l2a = 9a U) 12 = U) a 12 {R) = a 

With these assumptions we have 
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V S 
T t 
V t 

Vext 



a 



jJ2 $ + 2 * X! a "PaiU + ^PaiPa'ifafa' 
i a aa' 

'10 2 x -» 



by means of which the following equations of motions are obtained ( means transpose) 



/ Ml 5 
I map* mljv 



MwgZ 8 




tluj 2 In 



(Mg 
I N J \ 



(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 

(4.21) 



where X n is the identity matrix in n-dimensional space. 

If there are forces r l -" s ' ) between oscillators and the surface of the sphere (for instance thermal noise forces), then a 
new potential energy term must be added to the Lagrangian 



W ns = - ^ q t r. 



(ns) 



The equations of motion then become (cfr j?], eq. (25)]) 



/ MI 5 
I map* rnZfj 



—muj^ap 



X 5 ~ap 
In 



Mg 

r (ns) 



(4.22) 



(4.23) 



The coefficient a — ai2{R) can be calculated from (3.3) solving eigenfunction equations (2.4). One find that a — 
-2.886. 

Now the system of 5 + N coupled oscillators must be decoupled thus finding the eigenfrequencies. It is important to 
check if the frequency spread due to the decoupling is consistent with the assumption that only modes with n = 1 and 
I = 2 have been taken into account. Calling Co a these eigenfrequencies and f a the relaxation times of the eigenmodes, 
then the following consistency relation should hold 



u n i\ >> - t = min{f a ,T n i m } 

T 



n j£ 1 
I ^ 2 



As it will be seen, the above condition is fulfilled for high quality factors and rrii/M << 1. To perform such a program 
a particular arrangement of transducers must be chosen. It is what we are going to do in the next section. 



V. DODECAHEDRON ARRANGEMENT 



Now, the problem of the decoupling of the system made by the quadrupolar modes of the sphere interacting with 
N harmonic oscillators is to be solved. In order to do so the number and exact position of the mechanical oscillators 
on the surface must be known. This choice is of fundamental importance to keep isotropy in the sensitivity. Until 
now, two arrangements in particular are studied in view of an experimental realization. The first one |^,|7j is obtained 
placing six radial resonators along the directions of a dodecahedron (this arrangement was yet investigated [^| — p~0|| in 
the framework of the study of a network of resonant gravitational antennae). In the second one Q 5 transducers (four 
of which tangential) are placed in such a way that they could be excited by only one mode of the sphere; in this way 
decoupling is realized at once. In this paper we restrict ourselves to the first case. 

Let us consider therefore the dodecahedron arrangement. Five transducers are placed on the same parallel of the 
sphere. Their common azimuth angle is 9 — arcsin(2/\/5); their polar angles are given by 4>t — 2kTr/5 (k € {1, . . . , 5}). 
Last transducer is placed at a pole. In this case matrix p become 
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1 1 



C4 c 2 c 2 c 4 

-S 4 S 2 — »2 «4 

— 1 — C 2 — C4 — C 4 — Ca 

S 2 S4 — S4 — S 2 



\ 

* 






cos 



c 2 

,s 2 = sin 
a i 



'27^^ 

In \ 



-a\ —a\ —d\ —a\ —a\ 5ai / 
Such a matrix has the following important properties (cfr J5|,0]) 



5 

2\/3 



c 4 

S 4 



; (f) 



sin(f) 



(5.1) 



t ° T 

Z7T 



= 



(5.2) 



Decoupling system (4.23) means solving an eigenvalue problem. The eigenfunctions are the so called normal coor- 
dinates, while the eigenfrequencies are called frequencies of free vibration (see |l5[). In the case of this arrangement, 
normal coordinates are divided in three set: two quintuplet with eigenfrequencies lu + = ujq A + and c<j_ = ujq A_ and a 
singlet with frequency loq where 



A, 



1± 



3a^ 
2vr 



m 
M 



(5.3) 



The normal coordinates C p (p = 1 5 
C - Ce and C+ = (Ci) P € {7, . . 
appendix |A|) 



. . , 11) are therefore naturally divided in three groups: = (£ p ) p G {1, . . . , 5}, 
, . 11}. The equations of forced motion are (detailed calculations are found in 



(5.4) 



Therefore if g is given , on e can find the motion of coefficients f a and transducers qi. In particular because of the form 
of matrix A (see eq. (A9)) we have 



/ = 



2M ( A -<- 



A+Ch 



(5.5) 



'6m 



L (6xl) 



VI. DECOMPOSITION OF THE GRAVITATIONAL WAVE FORCE 



A gravitational wave, in the framework of a general metric theory of gravity, is described as a combination of six 
matrices. A possible choice is the following (see also [|l6||): 



V3 




(6.1) 



, , -1 \ / 1 \ / 1 ±i . 

#(o) = J_\ o -10 = ±i £;( ±2 ^ = ±t -1 ] (6.2) 

'002/ V 1 ±i / V000 



Matrix (6.1) accounts for the scalar part of the signal, while the others describe the quadrupolar components of the 
wave. They are relative to m — 0, m = ±1 and m — ±2 respectively. One may argue that, once the direction of 
propagation is known, a sphere with dodecahedron arrangement could detect, at least in principle, all six components. 
However, the scalar mode of the sphere has a principal resonant frequency u>io, which is a little more than twice the 
quadrupolar one uju (e.g. [Tq]). Current experimental schemes needs that, in order to measure the scalar mode of the 
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sphere, this must be coupled with a secondary oscillator having luiq as a resonant frequency. Therefore, being instead 
UJ12 the resonant frequency of each transducer, scalar part of the signal cannot be detected. Besides, propagation 
direction is in general unknown. One is led to the conclusion that with a dodecahedron arrangement it is possible 
to get the direction of propagation and three components of the five quadrupolar possible ones. In the framework of 
general relativity, a gravitational wave could only have two possible states of polarization, the E^ ±2 ^ ones. Therefore 
one can think that a sphere could be used to test general relativity. However, the discovery of PSR1913+16 O] and 
20 years of subsequent observations have put strong experimental limitations on other theory of gravitation different 
from general relativity. This last one passed all these new experimental tests with complete success flgf] showing that 
"the correct theory of gravity must make predictions that are asymptotically close to those of general relativity over a 
vast range of classical circumstances" [[l9| . Therefore another possibility is to consider general relativity as the correct 
metric theory of gravity and then to use the degree of freedom left as a veto against spurious signals. This is the 
strategy of the present work. In which follows we consider a gravitational wave in O' having the form 



E'{t) = E R {t)E 



(+2) 



E L {t)E 



(-2) 



Ei(t)E 



(0) 



(6.3) 



where El — E^ and the connection with the usual polarization amplitudes is E + (t) = En(t) + E^it) and E x (t) = 
i [En(t) — £^(i)), while Ei is the eventual longitudinal polarization. 

With this position it is possible to find, f or ev ery n, the vector g n (whose components are q n 2 a ) w hi ch e nters i n th e 
equations of motions of the normal modes (5.4) as a function of E'(t) by means of eqs. (4/7), ( 3.12 ), (2.S) and (|2.l|). 
Detailed calculations are found in appendix |B[ Here we give only the final result 



(gnhit) = -7= 



V2 



(9n)4{t) = 



- (3 - cos 2(5) cos 2H E + - sin 5 sin 2H Ey, 



1 + cos 25) cos 2H E t 



- (3 - cos 25) sin 2H E + + sin 5 cos 2H E y 




sin 25 cos H E< — cos 5 sin HEs, 



I - S1I 

sin 25 sin H E + + cos 5 cos HEy 



AT 

^- (1 + cos 26) sin2# Ei 

■ — sin 2(5 cos H E^j 

V3 \ 
: — sin 2(5 sin H E\ J 



/3 1 
^- (1 + cos2(5)£+ + - (1 -3cos2(5)£i 



(6.4) 



where j n is defined in appendix [B]. 

It has been therefore explicitly found the motion of mechanical transducers on the surface of the sphere as a function 
of direction and polarization amplitudes of a signal (by means of (|o~5|), (5.4) and ([O])), solving the direct problem of 
motion. 



VII. INVERSE PROBLEM 



Inverse problem lies in the determination of the incoming signal characteristics for a given response of the detection 
apparatus. In this case what is measured is the 6 component vector q(t), that is the motion of transducers with 
respect to the surface of the sphere. From this vector one should obtain C(£), desc ribi ng the motion of normal modes, 
and g, the generalized force for the normal coordinates. Finally, inverting system ( |6.4[ ) one should get those quantities 
which describe the signal, that is direction, polarizations E+ and Ey, and an eventual polarization with m = ±1 or 
to = 0. In which follows we consider only the possibility that the signal could have a longitudinal polarization Ei 
(see©). 

In Fourier space {f(o-i) = (2tt)^ 1 ^ 2 f(t) exp(~iujt) dt) using properties fl5.2| ) of matrix p and eq. ( |5.4p one gets 
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2tt M - 



2 



g{u) = -— ^ » 2 ' p*g( w ) (7.1) 



This last equation implies, together with (p.2| 



g*(w)g(o;) cx ff *(u;)flM = ^-{E\ + E 2 x + E?) (7.2) 



where last equality holds taking into account eqs. (p~7|), (3.12) and (Bl). Therefore the sum of the squares of the 



detector responses are independent on the direction. This means that this transducer arrangement gives isotropic 
sensitivity. The result could be also obtained calculating directly the sum of the squares of the components of g 



in ( |6.4| ). 

Inverse problem is then solved when system (|6.4|) is inverted. To this aim it is useful to define an amplitude vector 

a 

Pi 

a = y~g, a 2 = E\+ E 2 X + Ef, (7.3) 

7i 

whose length is the sum of the squares of the three polarization amplitudes of the signal. Performing this program 
(results are found in appendix |C|) one obtains a six-th degree algebraic equation for x = tan_ff/2 with coefficients 
depending on a. Then a relation for tan 5 as a function of H and a; finally three relations for E x , Ei as functions 
of H, S and a. 

As it will be seen better in the following, since for a gravitational signal Ei must vanish, only two of the six 
solutions for H can be accepted. In fact longitudinal components will be zero for only two values of H . These 
acceptable solutions differ only because they describe signals coming along opposite directions. In fact the difference 
between their hour angles is w. This property holds also for the discarded solutions: they are in pairs differing only 
for the opposite direction of propagation. 

Solution of the inverse problem lies therefore in the solution of the equation for H . This could be done numerically 
once vector a is known. To this purpose a code has been developed which solves the direct and the inverse problem. 
The obtained calculation precision on the reconstruction of the signal is very much higher than the errors introduced 
by thermal noise which we are going to consider in next section. 

VIII. THERMAL NOISE 

So far the theory dealt with ideal detectors, that is to say noiseless ones. However, as it is well known, real detectors 
are not at rest before the interaction with a signal but, because of the heat bath in which they are thermodynamically 
in equilibrium, they have a Brownian motion which give them a certain amount of energy. For this reason a detector 
(a part from the quantistic aspects of the problem, aspects which come out indeed when the temperatures are lower 
than those actually involved at the present) could only reveal signals yielding an amount of energy greater or at least 
comparable to the thermic one. Therefore, for any given temperature, there will be a lowest signal detectable. In this 
section we analyze how the behaviour of the system described until now is to be modified by the presence of thermal 
noise. 

To this aim, the problems that must be solved are at least two. The first one is to understand which is the noise that 
one could expect in the transducers on the surface of the sphere at a certain temperature. The second one concerns 
the knowledge of the decreasing in the precision with which not only the signal could be deconvolved (by means of 
the determination of the arrival direction and polarization amplitudes) but also with which could be used the veto on 
the longitudinal part of the signal. In other words a confidence level must be decided for which a certain value of Ei 
(which could not be exactly zero) is to be considered equal to zero or not. 

Let us begin to consider the first problem answering to the question: which noise is to be expected on the transducers 
at a given temperature? As a first estimate to be compared with the result of a more accurate calculation we may 
suppose, if the sphere mass is very much greater than that of the oscillators, that the Brownian noise cause a motion 
whose amplitude with respect the surface of the sphere is of the order of 



V muJ Q 

This result shall be valid as an order of magnitude, allowing settlements which, in general, depend upon the position 
of the oscillators on the surface of the sphere. We will see that, in the dodecahedron configuration this result is indeed 
correct being the same for every transducer. 







Strictly speaking, the equipartition energy theorem, that we used for giving the order of magnitude of the amplitudes 
Qi , is valid, in classical mechanics, only for every independent quadratic term in the energy and for every one of these 
the square-mean value is (1/2) fc^ T . Therefore, when we consider interacting oscillators we can apply it just to the 
normal coordinates. It must be observed, however, that the theory developed until now is only an approximation: in 
fact together with the lagrangian of the elastic body in interaction with the secondary oscillators we have not consider 
the dissipation function Jl5], §1-5]. As it is well known it is a quadratic function of the velocities T = (1/2) T^rfif 
and 2T is the rate of energy dissipation due to friction [115], §2.6]. The motion equations become [E5l §6.5] 







It is not in general possible to find normal modes for any arbitrary dissipation function. In other words it is not 
possible to decompose the system in non -interacting modes diagonalizing simultaneously the three matrices T^ , Vy 
and Tij. Therefore the equipartition theorem could not be applied even to the normal coordinates. If we make the 
assumption of high quality factor we can suppose that non-diagonal quadratic terms are very small with respect to the 
others so we can safely neglect them. In the framework of this approximation we can the refo re apply the equipartition 
theorem to the normal coordinates, whose Brownian motion will then be given by (see (5.4)) 



C- = -r n_(t) 



V2kT 
C = n(t) 



u> 
V2kT 



■2) 



n+(t) 



where n_(t), n(t) and n+(t) are eleven independent adimcnsional random functions (where "random" means that 
their average values vanish) whose root-mean-square deviat ions are equal to one. The relative coordinates q(t) are 
given as a function of the normal ones £(t) by the equation (5.5) which now could be written as 



V2kT 



— p* (n_(t)-n+(t)) + 
6m 



n(t) 
v6m 



L (6xl) 



1.3) 



Let us now define the following five functions 

v(t) 



V2 



(n-(*) -«+(*)) 



(8.4) 



For these functions (from now on < f > means the average of a given quantity /(i)) one obviously has that < v >= 
(that is to say they are random functions) and < v% >= 1. We can therefore write 



q(t) = - 

Let now be w(t) six functions defined as 



1 2kT 



3 m 



1 



/27rp*w(t) + l(6xi) 



(8.5) 



w(t) 



/ 2np t v(t) 



(8.6) 



As it is easily seen these are random functions and we have also that < wf >= 1. This follows directly from the fact 
that (p*p)ii = 5/(4tt). This property of the matrix p is of fundamental importance in order that all the transducers 
have the same noise amplitude. Introducing these new functions we obtain 



1 

w 



2kT 
3 m 



w(t) 



n(t) 1 



(6x1) 



•7) 



putting now 



z(t) 




-n{t) l (6xl) 



(8. 
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where once again < z >= and < zf >— 1, one gets 



q(t) = ±J—z(t) (8.9) 

u>o V m- 

We have therefore shown that the fluctuations due to the Brownian motion give rise to a random motion of the 



coordinates q(t) with root-mean-value just given by (8.1). It is to be remarked that this result depends on the 



particular arrangement of the secondary oscillators on the surface of the sphere. A priori one would in fact expect a 
position dependent result. 

We have therefore given the answer to the first question put by the presence of noise. The next section is devoted 
to the solution of the second problem envisaged here: the reconstruction precision. 



IX. NUMERICAL SOLUTION 



In this section we are interested in the answer to the problem of knowing the decrease in the signal reconstruction 
precision when brownian motion is to take into account. It is to be remarked that this precision decrease concerns not 
only the arrival direction and the transverse polarization amplitudes, but also the value of the longitudinal component 
which, in case of gravitational event, must be zero but, because of noise, could not be "mathematically" so. Therefore 
we have to precise the criterion with which this signal component is to be considered "physically" vanishing. Because 
of the complexity of the system we can not simply apply a kind of analysis based on the error propagation: the best 
solution is to perform numerical simulations. 

The strategy is to write a code which, once the signal characteristics are given, calculates the components of the 
vector a (that is to say the response of the device), gives them a gaussian noise with null mean value and unitary 
root-mean-square deviation and finally solves the inverse problem. When this is done six solutions are found which 
represents a possible signal. The four solutions which have the greater longitudinal polarization are discarded. With 
an iterative procedure over n attempts (n >> 1), one could extract two kinds of information (for every value of the 
signal-to-noise ratio): the first one is about the error on the reconstruction of the initial signal; the second one is 
about the typical value of the ratio between the longitudinal component Ei and the amplitude E = (E+ + E 2 ,) 1 ^ 1 ^ 2 ^ 
which one could expect from a real measure in the case of gravitational (that is to say transverse) signal. 

We have reached the conclusion that the signal could be reconstructed when the signal-to-noise ratio (SNR) is 
equal or greater than 10. If this is the case the relative error on the amplitude E and the absolute one on the 5 
angle is of the order of l/SNR. A great lack of precision appears on the H angle and on the polarization angle 
ip = arctan E x /E + (that is to say on the knowledge of the two polarization amplitude) when 5 assumes values closer 
to ±7r/2, as it was anyway well-founded to be expected. As far as the ratio Ei/E is concerned, it is distributed around 
with a typical spread which is of the order of the reciprocal of SNR. From this analysis one could conclude that, 
if the response of the detector lead to a reconstruction of the signal in which the six solutions are such that the ratio 
Ei/E is significantly greater than l/SNR, then this signal is certainly to be considered a spurious event and rejected 
as a possible gravitational event. 

In tables below we report some of the results and graphics concerning two simulations over 1000 attempts. In the 
first case the signal was supposed to come from a source placed at hour angle H = 1.0 rad and declination <5 = 0.7 rad. 
The amplitude E of the wave is 10 (in units of the SNR) , with the two polarization amplitudes given by E + = 8 and 
E x = 6. In the second case the source is placed at hour angle H = 1.0 rad and declination S = 1.3 rad. The signal 
amplitude is still given by E ~ 10, with E + = 6 and E x — 8. The results of the simulations are given in table ffl. The 
quantity Afl is the error on the solid angle; one has Afl = AH A5 cos S. In Figs. and | are shown the scattered 
plots of reconstruction position in the two cases. 

An analysis of the results shows immediately that the precision on H decrease when the source become nearer 
to the poles. However this have not substantially repercussions on the solid angle precision. Joined directly with 
the "spreading" of H is the uncertainty on the polarization amplitudes which increases when the signal direction 
approaches a polar one. 

A consequence of this analysis is that, once a reference frame O is chosen, the isotropy of the sphere with thermal 
noise concerns direction, total amplitude E and amplitude of the longitudinal component Ei of the signal. As far as 
the possibility of reconstruction of the polarization amplitudes E + and E x , it is strongly dependent on the declination 
angle. That is, sphere isotropy on the reconstruction of polarization angle tp = arctan(i? x /i? + ) is broken by thermal 
noise near polar singularity. 
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X. CONCLUSIONS 



An elastic sphere could represent, in a no longer far future, an important device in the context of gravitational 
astronomy. This lies in its peculiarity of total reconstruction of the signal with isotropic sensitivity and in the 
possibility of testing the transversality of the Riemann tensor. This last fact is very important because it allows to 
exert a veto against spurious signals. We think that the operation usefulness of a gravitational antenna does not 
depend only on its sensitivity or capability of signal reconstruction, but also on the confidence level with which an 
observed event could be associated to a gravitational wave (see ]8|-|To[ ) . The Lagrangian approach we have used, has 
been developed in the context of this strategy in order to write the correct equations of the coupled oscillators using 
the simpler formalism. In fact we think that a serious study should be done in order to find other configurations of 
the secondary oscillators for which also an eventual scalar amplitude of the signal could be measured and used as a 
veto. A first approach in this direction has been attempted with a different method by 

The analysis of the simulated output in the presence of thermal noise has shown that an amplitude signal-to-noise 
ratio is needed of almost 10, in order to exploit all the possibilities of a resonant sphere (cfr JToj] ) . 
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APPENDIX A: DECOUPLING 



In this appendix we decouple Eq. (1.21) in the dodecahedron arrangement. 
Let us set rj an 11-vector whose components are 

Vi = ft i € {1, . • . ,5} 

Vi+5 =qi i € {1, ...,6} 

With this position the Lagrangian of the system is written 



where 



Tij = M 



(1 + 3a 2 /x/27r)X 5 a^p 



MU 



is the kinetic energy matrix, while 



= Mw 2 



T 5 




is the potential energy matrix. In ( |A3| ) and (A4) we have put /i = m/M. Motion equation is 

Tirf + Vijrf = 

To decouple the motion, the eigenvalues and eigenvectors of the matrix 

(Vij - LU 2 Tij) = Mwl{vij ~ X 2 Uj) 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



must be found (see |15, Cap. 6]). In (A6) we set A = uj/lu . Putting equal to zero the determinant of matrix (A6) 
one obtains 



det(Vy - w 2 T y ) = (M^YV^l - A 2 ) 
The three distinct eigenvalues are therefore 



A ' 2 ( 1 + A j + 1 



(A7) 
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A ll2 ,3,4,5 = A_ = l--y— fi 

A fi = A = 1 



(A8) 



i 1 / 3a2 

^7,8,9,10,11 = A + = 1 + 2V "27^^ 



To look for eigenvectors of (A6) is t he same as t o find the matrix A which simultaneously diagonalizes the matrix of 
kinetic and potential energy T (|A3|) and V (A4). This matrix is found to be 



A 



1 

2M 



A_Z H 







(5x1) 



A + Z 5 



(A9) 



y /2Tr/3fi\-p t l (6x i)/V3/i -- v /27r/3//A+p t 



where we have introduced notation k^ Ly N ^ to indicate a matrix with dimensions L x N whose elements are all equal 
to k. A is also the transition matrix between the generalized coupled coordinates rji to the normal coordinates Q 



(A10) 



Let us now consider the forced system. Let F — (Fj) be the generalized forces corresponding to the generalized 
coordinates r] — (rjj). Then, generalized forces Q — (Qi) for the normal coordinates are 



Q = A l F 



If the external forces act only on the sphere modes then (see (4.21)) 



Ft = M 9l 
Ft = 



Therefore 



ie{i,. 
i e {6, . 



■ 5} 

■ 11} 



(All) 



(A12) 



If we put 



Q 



C = (Ci) 

C = Ce 

C+ = (CO 




(A13) 



ie{i,...,5} 



ie{7,...,ii} 



(A14) 



then eqs. (5.4) and (5.5) are immediately found. 



APPENDIX B: CALCULATION OF g 



In this appendix we show how to find the components of vector g n ( see Sec. VI). First the amplitudes R' n im'{t), 
defined in (2^) must be calculated. With the given gravitational signal ( |6.3|) one has that the only coefficients which 
does not vanish in O' are (cfr. with ||): 



^22(0 = 
^2-2(0 = ' 

-Rn2oM = 

where, using the same notation as [gj, we have 



32tt 
~15 
32^ 
HT 

167T 

~L5" 



R(a n2 + 30 n2 )E R (t) 



R(a n2 + 3/3 n2 )E L (t) 



R{a n2 + 3/3„ 2 ) Ei(t) 



(Bl) 
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* n2 = (MR)- 1 J a nl {r)pr 3 dr [3 n2 = M~ l J b nl {r)pr 2 dr 



(B2) 



Now G nm i(t) defined injpj), that is to say the coefhcients of the gravitational wave decomposition in O are to be 
evaluated. Because of (Bl) one has 



%l(t) - S12 (£>m2-ft'>22W + ^-2^2-2 (*) + I3 m0^n2o( < ) 



(H, 6) = exp [imH]d% m , (5 + tt/2) 



Following @J| one has 

where d|^, (/3) (/^ = ^ + 7r /2) are given by (see |14|]) 
d (2) (/?) = 



(2) 



/ d 22 (/3) -di 2 (/3) d 02 (/?) -d-i 2 (/?) d_ 22 (/?) \ 

dia(/3) diiGS) d 10 (/3) di_i(/3) -d_ ia (/9) 1 

d 02 (/3) -dio(/3) doo(/3) d 10 (/3) d 02 (/3) 

d_i 2 (/3) di-xO?) -di (/3) dn(/3) -d 12 (/3) 

Vd_ 22 (/?) d_i 2 (/3) d 02 (/3) d 12 (/3) d 22 (/3) / 



where 



(B3) 
(B4) 

(B5) 



d 22 {S + ir/2) = 3/8- 1/2 sin 5- 1/8 cos 25 
d^ 22 {5 + 7T/2) = 3/8 + 1/2 sin 5- 1/8 cos 25 

d 12 (5 + Tr/2) = -l/2coso" + l/4sin2<5 
d_ 12 (5 + 7r/2) = -l/2cos5- l/4sin2<5 

dnO + Tr/2) = -1/2 sin J - l/2cos25 
di_i(5 + 7r/2) = -1/2 sin 5 + l/2cos2o' 

do 2 (5 + vr/2) = v/3/32 (1 + cos 2(5) 

di (5 + 7r/2) = -V378sin2(5 

doo(<5 + 7r/ 2 ) = (1 - 3cos25)/4 

Setting 7„ = — -y/327r/15 R(a n 2 + 3/3 n2 ), the irreducible tensor of second rank G„, whose components are Gnim is 
written 

I A r,iH 

d - 12€ E t (t) 



( d 22 e aH \ 

«12 e 



E R (t) 



-i2e 
2 2e" 



-iH 
-i2H 



E L {t) 



J 



-d 



i 2 e 



( d 02 e aH \ 

A iH 
«10 e 



V2 



V d2 2e - i2ff / 



V d 02 e^ 2H J 



(B6) 



Applying to G n the operator f/*- 2 -' defined in (4.8) we get the real vector g n whose components are given by (6.4). 



APPENDIX C: INVERSION OF SYSTEM (6.4) 



After manipulation of system (6.4) with definition (7.3) we find the following equation for x — tanff/2 



C( ) + 2 C(i) x + C( 2 ) x 2 + 4 C( 3 ) x 3 - C( 2 ) x +2 C(i) a; 5 - c (0 ) £ 6 = 



(CI) 



where 



c (o) = a 2 a 3 ~ V^ai 
C(i) = a 2 3 - 2a 2 a 3 2 



a 2 0,4 — 03 04 + o 2 04 0,5 



+ V3 ai 03 04 + a 2 04 + 2V3 ai a 2 05 + 3 03 04 05 — 2 a 2 05 
C( 2 ) = -11 a 2 2 a 3 + 4a 3 3 + 

+ 7 V3 ai a 2 04 + 7 03 a4 2 — 8 \/3 a\ 03 05 — 15 a 2 0,4 05 — 8 03 as 2 
C( 3 ) = -a 2 3 + 4a 2 a 3 2 - 

— 3 V3 a\ 03 04 — 3 a 2 a 2 — 2 v3 ai a 2 05 — 5 03 04 05 + 2 a 2 a^ 2 
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Declination angle <5 is found from 



a 4 cos 2H - a 5 sin 2H 

tan <5= — — (C2) 

a,2 cos H — 0,3 sin ii 



while the three polarization amplitudes are given by 

E x = (a-4 cos 2H — 05 sin 2H) sin 5 + (02 cos H — 0,3 sin i?) cos 6 (C3) 

(#+ - VSEi^j = 2(a 3 cosH + a 2 sin if) sin 2(5 - (05 cos 2iJ + a 4 sin 2H - V3ai) cos 2<5 (C4) 

E+ + -=Ei ) = 05 cos 2fl" + o 4 sin 2F+-^= (C5) 



V3 V V3 
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FIG. 1. Simulation over 1000 attempts of reconstructed position of a source located at hour angle H = 1.0 rad and declination 
S — 0.7 rad. Amplitude signal to noise ratio is 10. The incoming gravitational wave is assumed linearly polarized with 
E x /E+ =3/4. 

FIG. 2. The same as Fig. | with H = 1.0 rad, 5 = 1.3 rad and E x /E + = 4/3. 
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TABLE I. Reconstruction over 1000 attempts of gravitational signals coming from sources placed at hour angle H and 
declination S. Afl is the angular resolution in steradians; E is the amplitude of the incoming waves in units of the SNR; E + 
and E x are the polarization amplitudes; Ei is the eventual longitudinal component of the signal, which is used as a veto against 
spurious events. 





First case 


Second case j 




Original data 


Reconstruction 


Original data 


Reconstruction 


8 (rad) 


0.7 


0.7 ±0.1 


1.3 


1.3 ±0.1 


H {rad) 


1.0 


1.0 ±0.15 


1.0 


1.0 ±0.5 


AQ (sterad) 




1 x 10^ 




1 x 10~ 2 


E (SNR) 


10 


10 ±1 


10 


10 ± 1 


E+ 


8 


8 ±1.5 


6 


4±5 


Ex 


6 


6±2 


8 


6.0 ±4.5 


Ei jE 





0.0 ±0.1 





0.0 ±0.1 
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